Critical behaviour of the XF-rotors model on regular and small world networks 
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We study the XY-rotors model on small networks whose number of links scales with the system 
size Nii n ks ~ N 1 , where 1 < 7 < 2. We first focus on regular one dimensional rings. For 7 < 1.5 
the model behaves like short-range one and no phase transition occurs. For 7 > 1.5, the system 
equilibrium properties are found to be identical to the mean field, which displays a second order 
phase transition at e c — 0.75. Moreover for 7 C = 1.5 we find that a non trivial state emerges, 
characterized by an infinite susceptibility. We then consider small world networks, using the Watts- 
Strogatz mechanism on the regular networks parametrized by 7. We first analyze the topology and 
find that the small world regime appears for rewiring probabilities which scale as psw oc 1/N 1 . Then 
considering the XY-rotors model on these networks, we find that a second order phase transition 
occurs at a critical energy e c which logarithmically depends on the topological parameters p and 7. 
We also define a critical probability pmf, corresponding to the probability beyond which the mean 
field is quantitatively recovered, and we analyze its dependence on 7. 
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I. INTRODUCTION 

Real-life networks are of finite size, loopy and display 
heavy correlations. This complexity represents a chal- 
lenge from several points of view: first it is computation- 
ally expensive when attempting to investigate the net- 
work topology and to simulate dynamical systems upon 
it; moreover it becomes rapidly intractable analytically 
and one is obliged to make assumptions in order to sim- 
plify the picture and perform calculations. If this effort 
is crucial to practically afford problems, it also embeds a 
deeper question: facing the network complexity and their 
omnipresence in real world, it is fundamental to make 
the distinction between the essential variables which are 
able to catch the topology main features and those de- 
tails which are unessential for a minimal though complete 
description. One of those very fruitful simplifications is 
sparseness, i.e. the networks considered have in general a 
few links per vertex while the network size tends to infin- 
ity. More precisely, a network is sparse if k/N — > when 
N — > 00, k being the average degree. This basic hypoth- 
esis leads to a crucial consequence: locally, the network 
can be approximated by a tree, which means the absence 
of finite loops, i.e. finite closed paths, among the vertices. 
Sparseness and the local tree-likeness proved essential to 
analytical studies of dynamical systems on networks: we 
cite, focusing just on small-world networks, studies on 
the Ising model [IIQ, percolation Q and, more recently, 
on the Kuramoto model (for a more complete overview, 
see [J]). Therefore, the advantage in terms of numerical 
computation is evident: in general, both numerical stud- 
ies investigating the network topology [g, Q and critical 
phenomena on networks [J Q, 0, B exploit the assump- 
tion of sparseness in its strongest form, taking the degree 
as constant. Nevertheless, increasing the links density, 
networks exist which are still sparse, fulfilling the afore- 
mentioned condition but they cannot no longer ensure 
the tree-likeness because of the heavy presence of loops. 
It could be argued hence that the links density could play 
a non negligible role both on the topological properties 



of those networks and on dynamical models defined upon 
them. Indeed it is the case of the XY-rotors model on 
regular one-dimensional chains: we show that the pas- 
sage between a sparse network in the sense of k = 0(1) 
and a dense one (k = O(N)) implies the emergence of a 
new metastable state for which the thermodynamic or- 
der parameter does not relax at equilibrium [9j . The links 
density hence triggers a non trivial effect on the thermo- 
dynamic behavior of the XY model, which by itself is 
known for possessing a rich phenomenology investigated 
in several numerical studies jl0NT5l | on 2 and 3-D lat- 
tices. In particular, we would like to recall, as an exam- 
ple among many others, that the two dimensional case 
with nearest neighbors coupling is characterized by the 
famous Berezinskii-Kosterlitz-Thouless phase transition 
[la . [l7fl , which implies the correlation function to switch 
from a power law decay at low temperatures to an ex- 
ponential one in the high temperatures regime. In the 
mean field limit, the XY model, called the Hamiltonian 
Mean Field (HFM) model in this case, displays as well a 
wide variety of behaviors this complexity being strongly 
entangled with the lack of additivity. Among its peculiar- 
ities we cite the presence of a second order phase transi- 
tion of the magnetization |18| and, even more noteworthy, 
the presence of non equilibrium quasi-stationary states of 
diverging duration in the thermodynamic limit |l9l— t22j| . 
More recently, those models have been challenged to face 
more complex network topologies: for instance, studies 
exist concerning the HMF model on random graphs [23J 
where, varying the links density, a second order phase 
transition of the global magnetization is recovered for 
every density value in the thermodynamic limit. Fur- 
thermore, studies of the XY model on small world net- 
works 0, [f| proved that this lattice topology supports as 
well complex thermodynamical responses of the model: a 
mean field transition of the order parameter is retrieved 
and its critical energy seems to depend on the network 
parameters. 

The present work inscribes itself on this line as we will 
focus, on first instance, on regular networks and then we 



will shuffle this regular topology with the introduction of 
a controlled amount of randomness. The first part of the 
paper being on regular networks, we detail the analyti- 
cal calculations presented in [9j showing that tuning the 
link density allows to pass from a short-range regime to 
a long-range one. The analytical approach is preceded in 
Sec. Illll bv the results of numerical simulations which are 
as well more extensively illustrated than in [9[. Further- 
more, we show that it exists between those two regimes 
a peculiar metastable state characterized by huge fluc- 
tuations of the order parameter. We then address, in 
the second part of the paper, small-world networks us- 
ing the Watts-Strogatz model |2J] aiming to shed light 
on the interplay between the link density and the injec- 
tion of randomness in the network. In his regime we 
first investigate, acting on the links density 7 and on 
the rewiring probability p, the crossover from the reg- 
ular lattice to the small-network topology. In Sec lIVBl 
we consider the dynamics of the Xy-rotors model on 
small-world networks and we show how the emergence 
of global coherence, via a mean field phase transition of 
the order parameter, strongly depends on the topological 
conditions fixed by p and 7. Furthermore, we discuss in 
the last part how this influence turns out to be quanti- 
tative, affecting the critical energy e c at which the phase 
transition occurs. 



II. THE XY-ROTORS MODEL 

The XY-rotors model describes a set of N spins in- 
teracting pairwise: each spin is fixed on the sites of a 
one dimensional ring and it is assigned with two canon- 
ically conjugated variables {6i,Pi}, Q\ G [— 7r; tt] being a 
rotation angle. The XY Hamiltonian reads [17ll25j: 
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are given by the set of Hamilton equations: 
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where Vi, represents the neighbors of rotor i. A global 
parameter, the magnetization is defined by 
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in order to have an insight on the macroscopic behavior: 
we expect finite values of M to indicate the emergence of 
a coherent inhomogeneous state, while a vanishing mag- 
netization signals the absence of long-range order. We 
first study the response of the total equilibrium magne- 
tization M to the change of the underlying network via 
the 7 parameter. Practically, for each 7, we perform 
simulations within the microcanonical ensemble, by di- 
rect numerical integration of Eqs. (jj) with the fifth order 
optimal symplectic integrator described in [2a|. The ini- 
tial conditions of angles and momenta are picked from 
a Gaussian distributions with identical variance (which 
corresponds to a low temperature setting) and, to check 
the numerical integration, we monitor the conservation 
of the two constants of motion preserved by the dynam- 
ics: the energy E = H and the total angular momentum 
P = J2iPii which we have set without loss of generality 
to P = 0. Finally the time step is At = 0.05 and we av- 
erage the thermodynamic quantities over time only when 
the system has reached the equilibrium. 



III. THERMODYNAMIC BEHAVIOR ON 
REGULAR LATTICES 

A. : Numerical Computation 



where ej ,■ is the matrix encoding the spins connections: 



1 if i, j are connected 
otherwise 



(2) 



We take J > 0, so that we are in the ferromagnetic case 
and in the following J = 1 as well as the lattice step. 
Finally the 1/fc factor in Eq. (TTJ) ensures that the energy 
is an extensive quantity, k is referred to as the degree 
and, to control the density of links in the network, we 
define it as: 
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Practically, we take the integer part of Eq. ([3]) since, once 
fixed 7 and N, k is in general non integer. The dynamics 



The regular network that we take into account is a 
one-dimensional chain of N spins (rotors) with periodic 
boundary conditions for which each spin is connected to 
its k nearest neighbors. By tuning the parameter 7, 
1 < 7 < 2 we act on the links density of the network. 
For 7 = 1 (k = 2) the spins are connected to their near- 
est neighbors, while for 7 = 2 (k = N — 1) the network 
is fully coupled. Heuristically, changing the value of 7 
corresponds to change the range of interaction of each 
spin. Then two limit behaviors naturally emerge from 
this approach: the first is 7 — >1 in which we expect the 
system to behave progressively like a one dimensional 
short range system with the existence of a continuous 
symmetry group, and so without any phase transition . 
On the other side, the 7 — > 2 limit leads to the mean field 
regime and we expect the HMF transition of the magne- 
tization to appear above a specific threshold of degree. 
We find this boundary value for 7 = 1.5 so that the two 
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FIG. 1: (a) Equilibrium magnetization versus energy density 
for 7 = 1.25 and different sizes. The error bars are of the 
size of the dots; (b) Residual magnetization for 7 = 1.25 at 
e = 0.1 versus the system size. 



aforementioned limits translate more precisely in two in- 
tervals 7 < 1.5 and 7 > 1.5. Practically, for each 7 value, 
we monitor the average magnetization M(N,e) (the bar 
indicates the temporal mean) for different sizes N and 
for every energy density e = H/N in the physical range. 
The temporal mean is computed on the second half on 
the simulations: we start with the Gaussian initial condi- 
tions described in Sec. Ulland we simulate the dynamics, 
calculating the magnetization at each time step. When 
the system reaches a stationary state for the magneti- 
zation, we take the temporal mean as the equilibrium 
value. 

We start our analysis with the 7 < 1.5 interval. The 
simulations are displayed in Figs.QJ,-b) and as mentioned 
the magnetization smoothly vanishes with the tempera- 
ture (Fig. [1^). We have to recall that for low tempera- 
tures, the magnetization can be non-zero as a finite size 
effect, so the results displayed should depend on the sys- 
tem size. This is confirmed in Fig. [1] where the trend 
for the magnetization to vanish with increasing the size 
is exhibited. To check with even larger sizes, we con- 
sider in Fig. [T|d the magnetization for a small energy 
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FIG. 2: (color online) Correlation function Cj for 7 = 1.25 
and N = 2 14 . 



density e = 0.1 and several sizes. The results clearly 
point out that the magnetization vanishes in the ther- 
modynamic limit. When looking at relaxation scales, 
we found that larger sizes took more time to relax to 
equilibrium. So typically in our simulations we take as 
final time tf = 20000 for sizes up to N = 2 16 and for 
N > 2 18 t f = 30000. 

Given these numerical results, we conclude that in the 
7 < 1.5 interval, the system is short-ranged and the 
Mermin- Wagner theorem applies imposing the order pa- 
rameter to vanish. Nevertheless, if long-range order is not 
possible, quasi long-range could still entail an infinite or- 
der phase transition of the correlation function, like in 
the two dimensional XY model with nearest neighbors 
interactions. We recall that this particular type of crit- 
ical phenomenon, first detected by Berezinskii, Koster- 
litz and Thouless |la|, is characterized by two different 
types of decay of the correlation function with distance: 
a power law or an exponential decay, respectively for low 
and high temperatures. In order to look for such possi- 
bility we computed the correlation function 
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[N]), 



for every e in the considered range. The results shown 
in Fig. [2] indicate that the decay behavior is also expo- 
nential for low temperatures, demonstrating the absence 
of the aforementioned phase transition. This could have 
been anticipated from the fact that finite size effects on 
the magnetization even though present were small, but 
possible tricky effects of the boundary conditions could 
come into play, so it was worth while checking. 

To summarize our result, we can conclude that, for 7 < 
1.5, the spin degree is still too low for the system to show 
long-range or quasi long-range. Interestingly, the short 
range behavior is still at play even for configurations like 
7 = 1.4 where each spin is under the influence of quite 
an important neighborhood since k ex N 0A in this case. 




FIG. 3: Equilibrium magnetization for N — 2 16 and different 
7. For 7 7^ 1.5 the error bars are of the size of the dots. 



Taking now in account the symmetric interval 7 > 1.5, 
the spins are connected enough to allow a coherent state 
to emerge: in Fig. [3] the magnetization undergoes a sec- 
ond order phase transition at s c = 0.75 which is well 
described by the HMF analytical curve. Again, around 
the delicate zone of the phase transition, finite size ef- 
fects induce a shift between the theoretical prediction of 
the HMF and the simulations, but they can be smoothed 
down increasing the size. We recall that this phenomenon 
is also present for the full coupling 7 = 2. As a conse- 
quence we then find that even with a degree remarkably 
inferior (e.g. for 7 = 1.6 ) than the full coupling con- 
dition, each spin possesses enough connections to trig- 
ger the global behavior of the system and give a finite 
magnetization (at low energies). Of course, in both the 
intervals 7 ^ 1.5, the equilibrium magnetization is still 
affected by fluctuations because of the finite size effects. 
To monitor these we measured the magnetization vari- 
ance a 2 = (M — M) 2 and we show in Fig.[5]that it scales 



with the system size like 



oc 1/N. 



(6) 



This scaling is the one expected for the equilibrium state 
thus confirming that the values in Figs. QJ,- H3 are repre- 
sentative of such state. 

Given the results presented in the previous discussions, 
a natural critical value appears, which characterizes the 
shift from the short range picture to the long range one: 
7 C = 1.5. We decided to investigate the system behavior 
at this critical threshold imposing 7 = j c . In fine, we ex- 
pect that the system will be in a peculiar state by itself 
which cannot be labeled as short or long ranged. Results 
are depicted in Fig. |3l We observe that for low energies, 
0.3 < e < 0.75, the averaged magnetization is finite even 
when increasing the size but it remains lower than the 
mean field value. The effect is clearer when we look at 
its temporal behavior. It is indeed totally different than 
in the other two regimes and the order parameter M 



shows large fluctuations which are orders of magnitude 
larger than for the other 7 regimes. We show in Fig. 3] a 
comparison a time series for the same energy and system 
size and different values of 7, namely 7 = 1.75 which dis- 
plays a finite magnetization with small fluctuations and 
the one 7 = 1.5, with large fluctuations. In order to con- 
trol the fact that these fluctuations are not an artifact 
of our initial conditions and that it is likely that the sys- 
tem does not relax on larger timescales than the previous 
configurations, we considered computation times up to a 
final time tf = 200000. Results are presented in Fig. 21 
where it appears that this regime with large fluctuations 
persists. We recall that for 7 gs 1-5 the simulation time 
was at most tf = 30000 and it was enough to reach a 
stationary state. Proceeding further, we notice that the 
amplitude of these fluctuations is not dependent on sys- 
tem size. We compare for instance N = 2 12 to N — 2 18 
in Fig. [l]c, d, and conclude that for the aforementioned 
energies there is no significant amplitude decrease with 
system size. More precisely, if we consider the variance 
a 2 as before, it appears that the scaling of the variance 
mentioned in Eq. ^ and coherent with in the 7 7^ 1.5 
regimes, is substituted by a flat behavior increasing N 
(see the results in Fig. [SJd). It is worth noticing that the 
influence the system size can be retrieved not in the fluc- 
tuations amplitude but in the typical fluctuation time 
scale. In Figs. HJ>d, it becomes obvious that fluctua- 
tions appear to slow down with the system size. This 
time-scale dependence on the size is reminiscent of out of 
equilibrium behavior in systems with long range interac- 
tions, namely the lifetime of the Quasi Stationary States 
(QSS) |20l . |27H29| and further investigations are ongoing 
to shed light on this effect and on its potential analogy 
with the HMF results. 

Heuristically, for 7 — j c it is like the if each spin does 
not possess enough connections to create a global order 
and establish the mean field but, nevertheless, the degree 
is sufficiently high (% = 1.5 corresponds to k = y/N) to 
avoid the vanishing of the order parameter in the ther- 
modynamic limit. The resulting behavior is reminiscent 
of a bistable regime oscillating between the M = con- 
figuration and the mean field value, which corresponds to 
a finite magnetization, we thus may expect some kind of 
intermittent behavior. In fine, the flatness of the variance 
suggests moreover that we observe a state with infinite 
susceptibility x considering its definition 



X ~ lim Na . 

JV->oo 



(7) 



To conclude our analysis in symmetry with the 7 ^ 1.5 
cases, we looked for a signature of this non trivial state 
in the correlation function but the fluctuations heavily 
affect it too so that it oscillates without showing a proper 
scaling. 
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FIG. 4: (color online) Time series for the magnetization with 
(a) N = 2 18 , e = 0.60; (b) N = 2 12 ; e = 0.44 ; (c) e = 0.44 
and (d) e = 0.52. 
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FIG. 5: (color online) Scaling of the magnetization variance 
(<r 2 ) with the size for 7 = 1.75 (a) and 7 = 1.5 (b). 



B. Analytical Calculation 

The numerical investigations illustrated point out that 
the degree triggers the shift from the pure one dimen- 
sional topology to the mean field frame. We now tackle 
this issue analytically aiming to retrieve the influence 
of the topology, encoded in the adjacency matrix ti t j 
(Eq. &§), in the thermodynamic properties. We thus 
compute the magnetization in the low energy regime and 
check if the correct behavior is recovered, namely a zero 
magnetization for 7 < 1.5 and a finite value for 7 > 1.5. 
At low energies we have a clear separation between the 
magnetization values, with M = and the mean field 
one, in which as e — >• M —} 1. In this limit, due to 
the ferromagnetic coupling it is natural to assume the 
difference 9i — 8j are small when a t j = 1 so that the con- 
nected spins are mostly aligned in order to minimize the 
free energy. We can hence develop the Hamiltonian at 
the leading order: 
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Ik 
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so that our system reduces to a collection of oscillators 
connected by e^-. We then choose to represent the spin 
field as a superposition of modes and following the recipe 
given in refs [13, [30| : 



^ = Eto 1 ^(i)cos(^+0O 
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(9) 



In Eq. ©, we sum over N modes so that the change of 
variables is linear and we observe that, given the period- 
ical border conditions, it just corresponds to perform a 
discrete Fourier transform. The amplitudes ai are, in our 
approach, the information carriers of the temporal behav- 
ior, hence the representation of the momenta pi is related 
to the one of the angles via the first Hamilton equation 
Pi = Qi. The phases <f>i are randomly distributed on the 
circle to ensure that the momenta pi are Gaussian dis- 
tributed in the limit N — > oo as theoretically predicted 
for the microcanonical ensemble. Following the approach 
described in [3(J, if we consider different sets of phases 
{4>i}m labeled by m we can interpret each set as a real- 
ization of the system, i.e. a trajectory in the phase space. 
Hence the process of averaging on random phases would 
correspond to ensemble averaging and lead to dynamic 
equations which, nevertheless, embed information about 
the thermodynamic state of the system, via the phase av- 
eraging. If we now inject Eq. © in the Hamiltonian © 
we obtain for the kinetic part K: 



N N \^ 2 / 4^ ' 



(10) 



where (...) stands for the average over random phases. In 
Eq. (JTDJ) we used the relation: 



(C0s(fci + 4>i)COs{kj + <f)j)) 



1 



For the potential, we have that the adjacency matrix e^j 
is a circulant one because of the definition of the regular 
network given in Sec|TIJ Hence we can diagonalize it and 
its spectrum {A;} is known: 



2^ 2nlj 

t£ cos (— ) 
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sin[(fc + l)lir/N] 
sm(ln/N) 



(11) 

where k is the spin degree of Eq. (J3]). To the leading 
order the potential will hence take the form: 
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In Eq. (TT2"]) we used the identity 
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which comes from the eigenvectors of a circulant matrix 
of size N are the columns of the unitary discrete Fourier 



transform ruatrix of the same size. We can then inject in 
Eq. (|12p the linear waves representation and average over 
the phases as we did for the kinetic part of the Hamilto- 
nian: 



01 

N 



W(i-x l) \eA = W(i-x l )c 



Having obtained the averaged Hamiltonian (H) = (K) + 
(V), we can deduce the averaged equation of motion, as 
anticipated, via the second Hamilton equation 



d fd(H) 



dt V dcti 



d(H) 
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and obtain 



on 



-(1 - Xi)ai = - 



wfai. 



(13) 



We have hence an equation for an harmonic oscillator 
whose frequency depends on the adjacency matrix spec- 
trum and, consequently, on the spin degree. We note 
that this approach is dependent from our low temper- 
atures approximation, but as mentioned we shall make 
use of this by the fact that depending on the value of 
7 we expect two clearly defined regimes of zero or finite 
magnetization. Our system is now completely encoded 
in terms of wave amplitudes {oti} and frequencies {u>i} 
which can be linked observing that, at equilibrium, we 
have the equipartition of the modes (pi's are Gaussian): 



r 4i>M 
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In order to compute M, we apply the same procedure, 
meaning that we average over the phases its expression 
given by Eq. ([5]) after we substitute the representation 
Eq. ©. We obtain [13 : 



(M) =IJ.7o(aj)(cos0o,siii0o), 



(14) 



where J is the zeroth order Bessel function and 9 is the 
average of the angles {9{\. This quantity is conserved 
because of the translational invariance, giving a constant 
total momentum P which is set at P = by our choice of 
initial conditions. As the final step to evaluate Eq. (fl"4|) . 
we recall that we are dealing with a low temperatures 
approximation so we can consider that the amplitudes 
af to be small at equilibrium and in the large system 
size limit [30|. This consideration allows to develop at 
leading order the product of the Bessel functions and, 
taking the logarithm of Eq. (Ti"4"|) , we finally obtain: 
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Eq. (|15p conjugates the thermodynamic information 
and the topological one because of the matrix spectrum. 




FIG. 6: (color online) Analytical magnetization (M) from 
Eq. (I15|l for T = 0.1 versus 7. Theory refers to the exact 
analytical solution of the HMF model. 
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FIG. 7: (color online) Path lengths starting from the blue 
vertex. 



work is fragmented in many parts [5j. It is noteworthy 
for the following to add that the rewiring injects mainly 
shortcuts of length of the order of the network size O (J\f) , 
so that a fine tuning of the interaction range by the means 
of the randomness p is not possible. 



From one side it actually realizes the purpose of matching 
these two levels of description; nevertheless the spectrum 
in Eq. (jlip carries now the system complexity, requir- 
ing Eq. (fT5)) to be evaluated numerically. In Fig. [5] we 
show, increasing the size, how this approximated expres- 
sion grasps the correct asymptotic behavior, giving the 
mean field value in the high 7 regime and vanishing for 
low 7. The transition becomes sharper at 7 C — 1.5 by 
increasing the size and gives hence confirmation of its 
critical signification as already pointed out by our nu- 
merical simulations of Sec lIHI 



IV. THE SMALL WORLD NETWORK MODEL 



In Sees. MI AjIITTbI we considered a regular chain as 
network topology and we illustrated how the degree 
drives the thermodynamic response of the XY-model on 
those lattices from the short-range regime to the long- 
range one. The natural following step to reorganize the 
topology is now to break the translational invariance of 
the regular chain previously considered and to introduce 
some randomness in how the spins are connected. For 
this purpose, we used the Watts-Strogatz model (W- 
S)j24| for small- world networks, which interpolates be- 
tween a regular network and a random one by the pro- 
gressive introduction of random lon g-ra nge connections. 
Following the algorithm devised in [2J|, each link is re- 
connected with probability p to a randomly chosen other 
vertex or is left untouched with probability 1 — p: long- 
range connections are hence introduced and the rewiring 
procedure injects disorder in the network since k, fixed 
by Eq. ^ at the beginning, is non-uniform afterward. 
The degree distribution decays exponentially since the 
rewiring is performed independently for every vertex |5j . 
Moreover, since ki sa (fc), a W-S network is not locally 
equivalent, even in the limit case of p = 1, to a random 
graph were eventually isolated vertex exist and the net- 



A. Network analysis 

The small- world regime embeds characteristics of both 
the regular lattice and the random network ones: the 
network keeps track of the initial configuration since, af- 
ter the rewiring, it still conserves a local neighborhood 
like a regular lattice; on the other hand the network ap- 
proaches, in the sense specified in Sec lIVl the random 
graph topology because of the shortcuts induced by the 
rewiring. In our context it hence emerges naturally the 
question of how the degree, which scales as k ~ iV 7_1 , 
could influence the scaling of topological quantities in 
competition with the rewiring probability p. For in- 
stance, a crucial passage in which the 7 parameter could 
play an important role is the crossover from the regular 
chain topology to the small-world regime. It is usually 
investigated by the scaling behavior of the average path 
length l(p, 7), defined as the average shortest distance 
between spins. This quantity has an algebraic increase 
I ~ N for a regular one-dimensional lattice with fixed de- 
gree fc, while for random networks it grows as I ~ log N. 
The passage between those two regimes is enhanced by 
the long-range connections which could allow the spins to 
behave coherently. Practically, since the network lacks a 
metric, the distance between two spins is calculated as 
the minimal number of edges to cross to go from one 
spin to the other, as shown in Fig. [7J To investigate the 
change between these two behaviors we perform numer- 
ical simulations, varying 7 and p: we use values for 7 
from 1.2 to 1.5 and p ranges from 10 -7 up to 10 -3 . N 
is fixed at 2 14 and we average over 10 network realiza- 
tions for each value of p. In Fig. [5£i we plot l(^i,p)l 7(7,0) 
versus 7. l(^,p) shows the known crossover behavior |24| 
but, considering the probability psw(N,p,~/) at which 
l("f,p) drops abruptly to the random network values, it 
appears evident that it is strongly dependent on 7. We 
have the following scaling for p$w (N, 7) [3| , using the 



degree definition in Eq. ([3]): 



1 



Psw 



N D kD 



(16) 



where D = 1 is the dimension of the initial regular lat- 
tice. In Fig. [8b we plot the estimation of psw{^,l) 
from the simulations versus 7 which effectively confirm 
the power law of Eq. (fT6|) . The degree is hence crucial to 
quantitatively determine the passage to the small-world 
regime; this dependence unveils its importance consider- 
ing that, on small- world networks, a correlation length 
can be defined [3| as 



i = i/{pkDf' D 



(17) 



and then psw in eq. (|16p is the probability of having 
£ = N . This is the key condition to achieve global co- 
herence and it clearly appears that the density of links, 
governed by the parameter 7, and the randomness in- 
jected by p concur in complexifying the network topol- 
ogy. In Sec. IIII1 we then move one step further dealing 
with the thermodynamics of theXF-rotors model on the 
small-world network and looking for the topological sig- 
nature of the 7 and p parameters in its properties. 



(a) 




B. 



Thermodynamic Behavior on Small World 
Networks 



In Sec. IIV Al we focused on the topological interplay 
of 7 and p parameters in establishing the small-world 
regime which, as explained, is noteworthy for its ambiva- 
lence, resembling both to a regular lattice and to a ran- 
dom graph. In this section we put the XY"-rotors model 
on a small-world network: the question we address now 
is to investigate the thermodynamic counterpart of the 
network complex topology. We focus the low 7 regime, 
i.e. 7 < 1.5. In this case we recall that the degree is still 
too low to induce long-range order by itself without the 
intervention of randomness and the network behaves like 
a one-dimensional chain. In the interval 7 > 1.5 the high 
degree already induces a mean field phase transition of 
the magnetization whose critical energy is e c = 0.75, as 
shown in Sec. Mil In this interval, even without the con- 
tribution of long-range connections, the network is con- 
nected enough to behave like a full-coupled one, which is 
the case of the Hamiltonian Mean Field model. On the 
other hand, in the case of random networks, it has been 
shown that the mean field phase transition appears for all 
7 > 1 [23]. We thus introduce progressively long-range 
connections with the rewiring probability p since, from 
Eq. (|17p . we expect to retrieve two regimes determined 
by 7 and p; £ < N in which long-range order is absent 
and £ > N where the order parameter displays a second 
order phase transition. In Figs. [9^,-c we set 7 = 1.25 
and, for each value of p, we consider several system sizes, 
above and below the threshold £(1.25, p,N) = N. The 
results displayed in Figs. [H] show the equilibrium mean 
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FIG. 8: (color online) (a) Average path lengths versus rewiring 
probability for different 7 values, (b) Power law scaling of 

psw(n)- 



value of the magnetization M versus the energy density 
e = E/N: for N = 2 12 , the probabilities p = 0.001 and 
0.005 (Figs.[5]a-b) are still too low to entail the crossover 
to the long-range regime and the system does not un- 
dergo a phase transition. On the other hand, the other 
two sizes considered N — 2 14 and 2 16 are in the £ > N 
regime and the mean field phase transition is recovered 
all the p taken in account. As explained, increasing the 
randomness decreases the small-world threshold; hence 
all the sizes show the phase transition of the magneti- 
zation for p = 0.01 (Fig. [Ht). Those results confirm 
the interpretation of £ in terms of a correlation length 
in the standard statistical physics sense: in Sec. IIV Al 
we showed that it signals the topological passage from 
regular to small-world network which identifies itself by 
a drop of the average path distance l(N,~/,p); equiva- 
lently in this short l(N, 7, p) regime the existence of long- 
range order is possible and, thus, we observe the second 
field phase transition of the thermodynamic order pa- 
rameter. Remarkably, the critical energy e c at which 
the transition occurs varies accordingly to the random- 
ness; we thus investigate this effect tuning 7 between 1.2 
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FIG. 9: (color online) Average magnetization M versus en- 
ergy density e = E/N for several system sizes and p — 
0.001 (a), 0.005 (6), 0.01 (c). 



and 1.5 and p from 10~ 7 to 10~ 3 . As explained before, 
it is worth focusing on the interval j < 1.5. In this 
case the shortcuts introduced by the rewiring process are 
crucial for the achievement of global coherence; while in 
the 7 > 1.5 we already know that phase transition with 
e c = Ehmf =0.75 occurs both on regular chains [9J and 
on random networks [2J|. In Fig. [TU1 we plot the critical 
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(b) 

FIG. 10: (color online) (a) Logarithmic dependence of the 
critical energy e c versus the rewiring probability p for different 
7 values, (b) Power law scaling of pmf versus 7. 



energy e c (p, 7) versus the rewiring probability p for sev- 
eral values of 7 and we observe that the phase boundary 
seems to be well described by the logarithmic form : 



log'(3(7> c ) 



(18) 



with C ~ 0.1. Eq. (|T5)) is coherent with the scaling pro- 
posed in [3, [H as far as the p dependence is concerned. 
Remarkably, in 0, Q , it was a result issued from Monte- 
Carlo simulations in the canonical ensemble while we 
work in the microcanonical frame. Moreover the afore- 
mentioned results of logarithmic scaling was found in the 
p — > regime, while where we are exploring regions with 
large values of p. We also have to insist on the fact that 
Eq. psp embeds an extra piece of information concern- 
ing the degree. Indeed, in our analysis the "quantitative" 
topological parameter 7 affects in his turn the critical 
energy e c through the function 3(7) , showing the non 
trivial role played by the links density in the thermody- 
namic behavior of the XF-rotors model. 

Another specific information can be retrieved from 
Fig. [TUa . There is a threshold beyond which a "satura- 
tion" process exists, to be more explicit for each value of 
7, we define a threshold probability pmf (7) for which the 
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critical energy is e c = 0.75, identical to values obtained 
in the mean field (7 = 2), or for the fully randomized 
networks [23|- For p > pmf(i) increasing the random- 
ness does not influence anymore the critical energy and 
in some way the resulting small- world network is, from 
the thermodynamic point of view, equivalent to a fully 
coupled graph. In Fig. 110b . we show how this probabil- 
ity threshold pmf(j) depends as a power law on the 7 
parameter. Note though that we expect that pmf ~> 
when 7 — > 1.5, because as we discussed before in the 
7 > 1.5 regime the system is already in the mean field 
state without any rewiring. Regarding Fig. 110b . we do 
not expect the results to be valid near 7 = 1.5. Indeed a 
precise estimation of pmf proves to be very delicate since 
it relies in its turn on the determination of the critical 
energy of the transition which is intrinsically a hard task. 
The simulations are performed with finite size systems. 
Therefore the measured pmf is influenced by these finite 
size effects and this depends also on on system size even 
we have actually pmf{i, N), we can then obtain a finite 
value also for 7 = 1.5. We also have to mention that 
we average also only on a finite number of network re- 
alization, which may also affect the results. Proceeding 
in our analysis, we recall that for the regular network 
a metastable state was found for j c = 1.5 in which the 
order parameter is affected by heavy fluctuations, sug- 
gesting that the system oscillates between low magne- 
tization values, proper of the 7 < 1.5 regime, and the 
mean field value of the 7 > 1.5 case. We can notice that, 
after the introduction of randomness, we do not observe 
this metastable state for 7 = 1.5 or any other value of 
7. In fact now a small (eventually vanishing) p is enough 
to generate a phase transition. It therefore exists an in- 
terplay between the "quantitative" parameter 7 and the 
"qualitative" parameter p; nevertheless those parameters, 
as anticipated in Sec. IIV1 are not equivalent when deal- 
ing with their influence on the thermodynamic behavior 
of the XY-model. This duality is so far not complete 
since it was not possible to retrieve the metastable state 
in the 7 < 1.5 regime acting exclusively on the p pa- 
rameter. In this sense, the randomness is "regularizing" 
the thermodynamic behavior: the rewired network sup- 
ports either the behavior of a regular lattice either, once 
the small- world regime is reached, gives rise to the phase 
transition of the magnetization. Summarizing, we can 
say that the noise created by the rewiring stabilize the 
passage between the two regimes and destroys the deli- 
cate metastable state which arose in the regular lattice. 



V. CONCLUSION 

In conclusion, we have studied the influence on the 
critical behavior of the XY-rotors model of two different 
network topologies, the regular lattice and small world 
network. In Sec. IIII1 we introduced the parameter 7 
which allows to tune the number of links from the lin- 
ear chain to the full coupling configuration. We identi- 



fied two main parameter regions: the first for 7 < 1.5 
in which the model has a one dimensional behavior and 
thus it does not display long-range or quasi long-range 
order as shown by numerical simulations. On the con- 
trary, in the second region (7 > 1.5), the spin degree 
is sufficiently high to lead to the emergence of a coher- 
ent state: we thus observe a mean field phase transition 
of the magnetization, identical to the one of the HMF 
model. More interestingly, we show numerical and an- 
alytical evidence of an unstable state at the threshold 
between the two regions, for 7 C = 1.5. In this peculiar 
state, the magnetization is affected by fluctuations which 
seem to be size independent and, furthermore, this state 
does not reach equilibrium on the timescales considered. 
We then calculated analytically an approximated expres- 
sion for the magnetization, obtained in the low tempera- 
tures regime, which demonstrates the topological critical 
nature of 7 C = 1.5. This expression retrieves correctly 
the two behaviors aforementioned and, since it contains 
the spectrum of the adjacency matrix, it points out that 
the topological origin of the three different phases shown 
by the simulations. We have then studied the role of the 
links density on the topology of small- world networks and 
its effect on the XY-rotors model dynamics. We have fo- 
cused, in Sec. IIV1 on the crossover to the small-world 
regime tuning the 7 parameter. We show by numerical 
simulations that psw has the scaling in Eq. (|16[) which 
is therefore consistent with J3|. Hence the links density, 
governed by 7, turns out to be crucial to enhance the 
crossover between the "large- world" regime and the small- 
world one cooperating with the rewiring probability p in 
the creation of long-range connections. We then investi- 
gated, in Sec. IIVB1 the thermodynamic response of the 
XY-rotors model to the variations of the network under- 
lying. We retrieved the emergence of a mean field tran- 
sition of the magnetization once p > psw- This latter 
condition implies the network to be in the £ > N case, 
using the definition in Eq. (|17|) . and thus confirms the 
interpretation of £ in terms of a correlation length from 
the statistical physics point of view. Moreover we found 
a logarithmic dependence of the critical energy £ c (p, 7) 
on p and 7 which lead to the scaling in Eq. (fTU)) . The in- 
terplay between the topological parameters in modifying 
e c saturates when e c = 0.75 which is the critical energy 
of the Hamiltonian Field model and hence we define a 
new threshold probability pmf which displays the power 
law scaling with 7 shown in Fig. [TUb . We find hence 
that a small (vanishing) amount of randomness regular- 
izes the 7 = 1.5 metastable state pointed out in Sec. IIIII 
and moreover it was not possible to recreate it in the 
7 < 1.5 interval just adding long-range connections with 
p. Therefore, as far as the thermodynamic behavior is 
concerned, we conclude that 7 and p are not equivalent 
when dealing with the transition to the mean field state; 
nevertheless, we anticipate here that a more refined cri- 
teria than randomness could be found in order to perturb 
the regular network in the low density regime (7 < 1.5) 
and enhance the creation of out-of-equilibrium effects like 



the 7 C = 1.5 metastable state. 
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